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Abstract 
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1 Introduction 

Hydrodynamic limit describes the time evolution (usually governed by a lim- 
iting PDE, called the hydrodynamic equation) of empirical density fields 
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in interacting particle systems (IPS). For usual models, such as the sim- 
ple exclusion process, the limiting PDE is a nonlinear diffusion equation or 
hyperbolic conservation law (see [20] and references therein). In this con- 
text, a random environment leads to homogeneization-like effects, where an 
effective diffusion matrix or flux function is expected to capture the effect 
of inhomogeneity. Hydrodynamic limit in random environment has been 
widely addressed and robust methods have been developed in the diffusive 
case (HH H2l HH HH [2H [261 121])- 

In the hyperbolic setting, due to non-existence of strong solutions and non- 
uniqueness of weak solutions, the key issue is to establish convergence to 
the so-called entropy solution (see e.g. [30]) of the Cauchy problem. The 
first such result without restrictive assumptions is due to [28] for spatially 
homogeneous attractive systems with product invariant measures. In ran- 
dom environment, the few available results depend on particular features of 
the investigated models. In [6], the authors consider the asymmetric zero- 
range process with site disorder on Z d , extending a model introduced in |10j . 
They prove a quenched hydrodynamic limit given by a hyperbolic conser- 
vation law with an effective homogeneized flux function. To this end, they 
use in particular the existence of explicit product invariant measures for the 
disordered zero-range process below some critical value of the disorder pa- 
rameter. In [22], extension to the supercritical case is carried out in the 
totally asymmetric case with constant jump rate. In [2S], under a strong 
mixing assumption, the author establishes a quenched hydrodynamic limit 
for the totally asymmetric nearest-neighbor i^-exclusion process on Z with 
site disorder, for which explicit invariant measures are not known. The last 
two results rely on a microscopic version of the Lax-Hopf formula. However, 
the simple exclusion process beyond the totally asymmetric nearest-neighbor 
case, or more complex models with state-dependent jump rates, remain out- 
side the scope of the above approaches. 

In this paper, we prove quenched hydrodynamics for attractive particle sys- 
tems in random environment on Z with a bounded number of particles per 
site. Our method is quite robust with respect to the model and disorder. We 
only require the environment to be ergodic. Besides, we are not restricted 
to site or bond disorder. However, for simplicity we treat in detail the mis- 
anthropes' process with site disorder, and explain in the last section how 
our method applies to various other models. An essential difficulty for the 
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disordered system is the simultaneous loss of translation invariance and lack 
of knowledge of explicit invariant measures. Note that even if the system 
without disorder has explicit invariant measures, the disordered system in 
general does not, with the above exception of the zero-range process. In par- 
ticular, one does not have an effective characterization theorem for invariant 
measures of the quenched process. Our strategy is to prove hydrodynamic 
limit for a joint disorder-particle process which, unlike the quenched process, 
is translation invariant. The idea is that hydrodynamic limit for the joint 
process should imply quenched hydrodynamic limit. This is false for limits 
in the usual (weak) sense, but becomes true if a strong hydrodynamic limit 
is proved for the joint process. We are able to do it by characterizing the 
extremal invariant and translation invariant measures of the joint process, 
and by adapting the tools developed in [5] . 

The paper is organized as follows. In Section [2], we define the model and state 
our main result. Section[3]is devoted to the study of the joint disorder-particle 
process and characterization of its invariant measures. The hydrodynamic 
limit is proved in Section |H Finally, in Section we consider models other 
than the misanthropes' process: We detail generalizations of misanthropes 
and fc-step exclusion processes, as well as a traffic model. 

2 Notation and results 

Throughout this paper N = {1,2,...} will denote the set of natural num- 
bers, and Z + = {0, 1, 2, ...} the set of non- negative integers. The integer part 
[x\ G Z of x G R is uniquely defined by |_^J < x < [x\ + 1. We consider 
particle configurations on Z with at most K particles per site, K G N. Thus 
the state space of the process is X = {0, 1, • • ■ , K} z , which we endow with 
the product topology, that makes X a compact metrisable space, with the 
product (partial) order. 

The set A of environments is a compact metric space endowed with its Borel 
cx-field. A function / defined on A x X (resp. g on A x X 2 ) is called local 
if there is a finite subset A of Z such that /(a, 77) depends only on a and 
(rj(x),x G A) (resp. g(a,r],£) depends only on a and (r](x),£(x),x G A)). 
We denote by r x either the spatial translation operator on the real line for 
x G R, defined by r x y = x + y, or its restriction to x G Z. By extension, if / is 
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a function defined on Z (resp. R), we set r x / = / or x for 16Z (resp. R). In 
the sequel this will be applied to particle configurations r\ G X, disorder con- 
figurations a G A, or joint disorder-particle configurations (a, 77) G A x X. 
In the latter case, unless mentioned explicitely, t x applies simultaneously to 
both components. 

If r x acts on some set and \i is a measure on this set, t x \i — ii o r" 1 . We 
let A^ + (R) denote the set of nonnegative measures on R equipped with the 
metrizable topology of vague convergence, defined by convergence on con- 
tinuous test functions with compact support. The set of probability mea- 
sures on X is denoted by 'P(X). If 77 is an X- valued random variable and 
v G P(X), we write rj ~ v to specify that 77 has distribution v. Similarly, for 
a G A, Q G V(A), a ~ Q means that a has distribution Q. 

A sequence (i/ n , n G N) of probability measures on X converges weakly to 
some v G V(X), if and only if lim^oo J f dv n = J f dv for every continuous 
function / on X. The topology of weak convergence is metrizable and makes 
"P(X) compact. A partial stochastic order is defined on "P(X); namely, for 
/ii,/i2 G V(X), we write /ii < \ii if the following equivalent conditions hold 
(see e.g. [2U EI]): (i) For every non- decreasing nonnegative function / on 
X, f f d\i\ < f f d\i2- (ii) There exists a coupling measure /Z on X x X with 
marginals \i\ and fj,2, such that ~fi{(r},£) : 77 < ^} = 1. 

In the following model, we fix a constant c > and define A = [c, l/c] z 
to be the set of environments (or disorders) a — (ce(x) : xfZ) such that 

Va; G Z, c < a(x) < c" 1 (1) 

For each realization a G A of the disorder, the quenched process (r)t)t>o is a 
Feller process on X with generator given by, for any local function / on X, 

L a f{v) = a ( x Wy ~ x )Kv(x)Mv)) If (v x > y ) - f(v)\ (2) 

where rj x ' y denotes the new state after a particle has jumped from x to y (that 
is rf ,v {x) = r)(x) — 1, r] x,y (y) = r)(y) + 1, r] x ' y (z) = r)(z) otherwise), the parti- 
cles' jump kernel p is a probability distribution on Z, and b : Z + x Z + — > R + 
is the jump rate. We assume that p and b satisfy: 
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(Al) Irreducibility: For every z G Z, Ylnen\P* n ( z ) + P* n (~ z )} > 0; where 

*n denotes n-th convolution power; 

(A2) finite mean: J2 z ^z \ z \p( z ) < +°°; 

(A3) fsT-exclusion rule: 6(0, .) = 0, 6(., K) — 0; 

fA^) non-degeneracy: 6(1, K — 1) > 0; 

(A5j attractiveness: 6 is nondecreasing (nonincreasing) in its first (second) 
argument. 

For the graphical construction of the system given by ([2]) (see [5] and ref- 
erences therein), let us introduce a general framework that applies to a 
larger class of models (see Section below). Given a measurable space 
(V, J-y,m), for m a finite nonnegative measure, we consider the probabil- 
ity space (Q, J 7 , P) of locally finite point measures u(dt, dx, dv) on IR + x Z x 
V, where T is generated by the mappings u u>(S) for Borel sets 5* of 
R + x Z x V, and P makes u a Poisson process with intensity 



denoting by A either the Lebesgue or the counting measure. We write E for 
expectation with respect to P. For the particular model ([2]) we take 

V := Z x [0, 1], v = (z,u) e V, m(dv) = c- 1 \\b\\ oo p(dz)X [07l] (du) (3) 

Thanks to assumption (A2), for P-a.e. u, there exists a unique mapping 



satisfying: (a) t \-t r] t (a,r] ,u) is right-continuous; (b) r] (a,r] ,u) = r] ] (c) 
the particle configuration is updated at points (t, x,v) G u (and only at such 
points; by (t, x,v) G u we mean u{(t,x,v)} = 1) according to the rule 



M(dt, dx, dv) = Ak+ 



(dt) Xz(dx)m(dv) 



(a, r)o, t) G A x X x IR+ ^ r) t = r) t (a, r) Q , u) G X 



(4) 



rjt{a, 770, u) = T a ' x ' v r] t - (a, r) Q , u) 



(5) 



where, for v = (z, u) G V, T a,x,v is defined by 




otherwise 



b(rj(x),r](x + z)) 



(6) 



Notice the shift commutation property 



q-r x a,y,v _ <j-i 
' 1 x — ' x I 



(7) 
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where r x on the right-hand side acts only on 77. By assumption (A5), 

j-a,x,v : x — ^ X is nondecreasing (8) 

Hence, 

(a, 7/0, £) H- r)t(a,r)o,uj) is nondecreasing w.r.t. 7/0 (9) 
For every a G A, under P, (?7t(a, 770, w))t>o is a Feller process with generator 

Laf(v) = E / ^ V"™*) ~ /(»?)] ( 10 ) 

With (EJ), (TTU|) reduces to (12]). Thus for any t G R + and continuous function 
/ on X, JE[f(r] t (a, 7/ , lo))} = S a (t)f(rj ), where S a denotes the semigroup 
generated by L a . From fl9]), for /ii, ^2 £ ^(X), 

A*i < A*2 Vi G R + , /ii5 a (t) < /^(t) (11) 

Property (fTT|) is usually called attractiveness. Condition (JS} implies the 
stronger complete monotonicity property ([131 EDj that is, existence of a 
monotone Markov coupling for an arbitrary number of processes with gener- 
ator (j2J), see (T57j) below; we also say that the process is strongly attractive. 

Let iV G N be the scaling parameter for the hydrodynamic limit, that is, 
the inverse of the macroscopic distance between two consecutive sites. The 
empirical measure of a configuration 7/ viewed on scale N is given by 

vr iV (7 ? )(rfx) = N- 1 Y^v(y)S y /N(dx) e M + (R) 

where, for x G R, 5 X denotes the Dirac measure at x. 
Our main result is 

Theorem 2.1 Assume p(.) has finite third moment. Let Q be an ergodic 
probability distribution on A. Then there exists a Lip schitz- continuous func- 
tion on [0,K] defined below (depending only onp(.), &(.,.) and Q) such 
that the following holds. Let (t)q, N eN) be a sequence of X.-valued random 
variables on a probability space (fi , J-o, Po) such that 

lim n N (riF)(dx) = u L)dx P -a.s. (12) 

N— >oo 
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for some measurable [0, K]-valued profile Uq(.). Then for Q-a.e. a G A, the 
Po<S>P-a.s. convergence 

lim n N (r] Nt (a, t)q (u ),u)))(dx) = u(.,t)dx 

N— >oo 

holds uniformly on all bounded time intervals, where (x, t) h- > t) denotes 
the unique entropy solution with initial condition uq to the conservation law 

d t u + d x [G Q (u)] =0 (13) 

We refer the reader (for instance) to |30j for the definition of entropy solu- 
tions. To define the macroscopic flux G®, let the microscopic flux through 
site be 

j(oi,v) = 3 + ( a >v) ~ j~(<*,v) ( 14 ) 

j/.zfcX: y<l)<y \ z 
y,z&L: y+z<0<y 

We will show in ^orollary 13.11 below that there exists a closed subset UP of 
[0, if], a subset A Q of A with Q-probability 1 (both depending also on p(.) 
and &(., .)), and a family of probability measures [y^ p : a G A Q , p G 7Z Q ) on 
X, such that, for every p G 7Z Q : 
(Bl) For every a G A^, v® ,p is an invariant measure for L a . 
(B2) For every a G A Q , ^ ,p -a.s., 

lim(2/ + l)" 1 V r](x)=p 

l-^oo * — ' 

xSZ: |a;|<Z 

(7?5) The quantity 

G?(P):= J J(a,v)^' p (dv) (15) 

does not depend on a G A Q . Hence we define G Q (p) as f|T5|) for p G 7?. Q and 
extend it by linear interpolation on the complement of TZ9, which is a finite 
or countably infinite union of disjoint open intervals. 

The function G® is Lipchitz continuous (see Remark 13.31 below), which is 
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the minimum regularity required for the classical theory of entropy solutions. 
We cannot say more about in general, because the measures v^' p are not 
explicit. This is true even in the spatially homogeneous case a(x) = 1, unless 
b satisfies additional algebraic relations introduced in [Sj. In the absence of 
disorder, for the exclusion process and a few simple models satisfying these 
relations (see for instance [3] Section 4]), we have an explicit flux function. 
Nevertheless, invariant measures are no longer computable when introducing 
disorder, so that the effect of the latter on the flux function is difficult to 
evaluate. However, in the special case b(n,m) = l{ n>Q }l{ m< K}, p(l) = 1, 
is shown to be concave in [29], as a consequence of the variational approach 
used to derive hydrodynamic limit. But this approach does not apply to the 
models we consider in the present paper. 



3 The disorder-particle process 

In this section we study invariant measures for the markovian joint process 
(®uVt)t>o on A x X with generator given by, for any local function / on 
A x X, 

Lf(a, V)=Y, f[f (a, 7™77) - f(a, rj)] m(dv) (16) 
that is, for the particular model 

L f(u, v)= Yl a ( x )p(y ~ x )Kv(x), v(y)) If («, v x ' y ) - /(<*, V)] ( 17 ) 

x,y(^E 

We denote by (S(t),t G M + ) the semigroup generated by L. Given a = a, 
this dynamics simply means that a t = a for all t > 0, while (rit)t>o is a 
Markov process with generator L Q given by ([2]). Note that L is translation 
invariant, that is 

r x L = Lt x (18) 

where r x acts jointly on (a,rj). This is equivalent to a shift commutation 
property for the quenched dynamics: 

where, since L a is a Markov generator on X, the first t x on the r.h.s. acts 
only on 77. We need to introduce a conditional stochastic order. For the 
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sequel, we define the set O = + U O-, where 



+ = {K^jGAxX 2 :,,^} 

O- = {M,()6AxX 2 :(<,} (20) 



Lemma 3.1 For two probability measures /i 1 = n l (da,dr]), fi 2 = fi 2 (da,dr)) 
on A x X, the following properties (denoted by fj, 1 <C fi 2 ) are equivalent: (i) 
For every bounded measurable local function f on A x X, such that f(a, .) is 
nondecreasing for all a & A, we have J f dfi 1 < J f dfi 2 . (ii) The measures 
fi 1 and fi 2 have a common a -marginal Q, and fi 1 (di]\a) < {i 2 (drj\a) forQ-a.e. 
a G A. (Hi) There exists a coupling measure Jt(da,dr],d^) supported on + 
under which (a,rj) ~ fi 1 and (a, £) ~ /i 2 . 



Proof of lemma HOI (ii)^(i) follows from conditioning. For (i)^(ii), con- 
sider f(a,rj) = g(a)h(r]), where g is a nonnegative measurable function on 
A and h a nondecreasing local function on X. Specializing to h = 1, using 
both / and — / in (i), we obtain 

g(a)fi 1 (da, drj) = / g(a)fi 2 (da, drj) 



Thus fi 1 and /j, 2 have a common a- marginal Q. Now with a general h, by 
conditioning, we have 



h{rf) fi 1 (dr]\a) j Q(da) < j g(a) I J h(i])fi 2 (dr]\a) J Q(da) 
Thus, for any nondecreasing local function h on X, 

h{rj) fj, 1 (drj\a) < / h(rf) [x 2 (drj\a) 



holds Q(da)-a.e. Since the set of nondecreasing local functions on X has a 
countable dense subset (w.r.t. uniform convergence), we can exchange "for 
any h" and "Q-a.e." In other words, ii l (dr)\a) < ii 2 {dr]\a) for Q-a.e. a e A. 

For (ii)^(iii), by Strassen's theorem ([31]), for Q-a.e. a G A, there exists 
a coupling measure ~jl a (dri, d£) on X 2 under which i] ~ fj}(.\ct), £ ~ /i 2 (.|a), 
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and r] < £ a.s. Then p(da,drj,d^) := f A [5p(da)p a (dr], d^)}Q{df3) yields the 
desired coupling. (iii)^(i) is straightforward. □ 

We now state the main result of this section. Let Z^, S and S A denote 
the sets of probability measures that are respectively invariant for L, shift- 
invariant on A x X and shift-invariant on A. 

Proposition 3.1 For every Q G <Sg , there exists a closed subset 7lQ of 
[0, K] containing and K , such that 

(l L nS) e = {^,Q6 5 e A , P G^} 

where index e denotes the set of extremal elements, and (v® ,p : p G 7l9) is a 
family of shift-invariant measures on A x X ; weakly continuous with respect 
to p, such that 



r](0)iy Q ' p (da,d7]) = p (21) 
lim(2/ + l)" 1 V rj(x) = p, u Q ' p - a.s. (22) 



l— ¥00 

x£7,:\x\<l 

p< p > =>. U Q ' P < v Q ' p ' (23) 

For p = G Vfi (resp. p = K G we get the invariant distribution 

5q Z (resp. <5f^), the deterministic distribution of the configuration with no 
particles (resp. with maximum number of particles K everywhere). 

Remark 3.1 The set 7l9 and measures z/^' p also depend on p(.) and b(., .), 
but we did not reflect this in the notation because only Q varies in Proposition 

m\ 

Corollary 3.1 (i) The family of probability measures v&{-) '■= v® ,p {.\ot) 
on X satisfies properties (B1)-(B3) on page\^ (ii) for p G 1Z Q , G Q (p) = 
J j(a, r\)vQ' p {doi, dr\). 

Remark 3.2 By (ii) of Corollary and shift-invariance of v^ ,p {da,dri), 
G Q (p) = [ j(a,rj)u^ p (da,dr}) = [ J(a,rj)u^ p (da,dr]) (24) 
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for every p G 1Z® , where 

j(a, rj) : = a(0) ^ zp(z)b( V (0), r)(z)) (25) 

Thus one can alternatively take J(a,r]) as a microscopic flux function (we 
refer to p. 1347] for an analogous remark in the non- disordered setting). 



Proof of corollary \3.1\ Properties (Bl) and (B2) follow from Proposition 13. II 
by conditioning (here and after, we proceed as in the proof of Lemma 13. ip . 
By translation invariance of v®' p and conditioning we have, for Q-a.e. a G A, 



T4%* = p?£ (26) 

where t x on the l.h.s. acts on X. For property (B3) the result will follow 
from ergodicity of Q once we show that, for every p G 1Z Q , G®(p) = G® ia (p) 
holds Q-a.s. To this end we note that, as a result of f fl9|) . 

Laivi 1 )} =j(<x,v) -j(n<x,nr)) 

Taking expectation w.r.t. invariant measure v® ,p , and using (1261) . we obtain 
Gg(p) = f j(a, V )u^(d V ) = [ jinocnv^idr)) 

□ 

To prove Proposition 13. 1[ we need some definitions and lemmas. For ev- 
ery a G A, we denote by L a the coupled generator on X 2 given by 



(27) 



for any local function / on X 2 . For the particular model (jBJ), this is equivalent 
to the "basic coupling" of L a defined in [8], namely L a = Ylxyez-x^y^a > 
with 7% v f(ri,€) given by 

a(x)p(y - x) [b( V (x),rj(y)) A 6(£(x), £(y))] [f(v x ' v , C' y ) ~ f(v, 0] 
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+ a(x)p(y - x) [b( V (x),r](y)) ~ &(£(*), ttvW lf(v x ' v , - f(V, 0] 
+ a(x)p(y - x) [b{£(x),Z{y)) - b( V (x), V (y))}+ [f( V , C' y ) ~ f(v, 0] (28) 



If (r} t ,£t) is a Markov process with generator L a , and r] < £ > then rjt < it 
a.s. for every t > 0. We indicate this by saying that L Q is a monotone 
coupling of L a . We denote by L the coupled generator for the joint process 
(a t , r) t , 6)t>o on A x X 2 defined by 



for any local function / on A x X 2 . Given = a, this means that a t = 
a for all t > 0, while (?7 t ,£t) t >o is a Markov process with generator L a . 
Let S(t) denote the semigroup generated by L. We denote by S the set of 
probability measures on A x X 2 that are invariant by space shift r x (a, T], £) = 
(r x a, t x T], t x £). In the following, if V(da, drj, d£) is a probability measure on 
A x X 2 , Vi, T>2 and z7 3 (resp. z7 12 and z7 13 ) denote distributions of a, 77 and £ 
(resp. (a, 77) and («,£)) under 77. 

Lemma 3.2 Let G (X £ fl S) e with a common a-marginal Q. Then 

there exists V G (X^ D5) e sttc/i /j/iai 77 12 = /i' and 77 13 = fj," . 

Proof of lemma WTB. Let .M (//, //') denote the set of probability measures 
V G 1^(1 5 with z7 12 = // and 77 13 = //'. We show that .M (//,//") is a 
nonempty set. Set lP{da,dr],d^) := Q(da)[fi'(dr]\a) <g> ji"(d£,\a)). Then zv^ = 
//, 17O3 = /i" and z7° g 5. Let 



The set {z7*,t > 0} is relatively compact because v\ = Q is independent of t 
and, for i G {2, 3}, v\ < 5^ z . Let V°° be any subsequential weak limit of V t 
as t — > 00. Then 77°° retains the above properties of 77°, and 17°° G X^, thus 
77°° G Ai(n',fi"). Let /7 be an extremal element of the compact convex set 
-M(//, //'). We now prove that 77 G (Zj; fl S) £ . Assume there exist A G (0, 1) 
and probability measures V 1 , V r on A x X 2 , such that 



Lf(a,r],0 = (L a f(a,.))(r),0 



(29) 




v = \v + (1 — \y 



(30) 
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with v % G Tj; PI S for z G {7, r}. Since z/ G M.{p! ', //') , the projections of (1311 
on (a, 77) and (a, £) yield 

// = A^ 12 + (1-AK 2 (31) 
// = + (32) 

For i G {/, r}, V 1 G X^ D «S implies u[j G X^ D S for j G {2, 3}. Since //, //' 
belong to (X L n «S) e , z7* 12 = //, = // by (}3"T]) - (l32"l) . that is, V 1 G AT(//, //')■ 
Since V is extremal in Ai(fJ,', //'), ( 1301 yields V 1 = V r = V. □ 

Lemma 3.3 Lei v be a stationary distribution for some Markov transition 
semigroup, and (X t ) t > a Markov process associated to this semigroup with 
initial distribution v. Assume A is a subset of E such that, for every t > 0, 
l A (X t ) > Ia(-^o) almost surely. Then v A (dx) = v{dx\x G A) and v A c{dx) = 
v(dx\x G A c ) are stationary for the considered semigroup. 



Proof of lemma \3.3[ Since (X t ) t > is stationary, we have E[l y i(X t )] = 
E[lyt(Xo)], thus 1a(X q ) = l A (X t ) and l A o(X t ) = l A c(X ) almost surely. 
Stationarity of u A amounts to TB[f(X t )\X G A] = JE,[f(X )\X G A] for 
every bounded /. We conclude with 

mf(xMY ^ M nf(Xt)lA(X )} E\f(X t )l A (X t )] 

E[/(X )1 A (X )] wIHxMXgA] 

□ 

Lemma 3.4 Let V G (l E n5) e . ThenV (0 + ) andu(0^) belong to {0,1}. 



Proof of lemma [5^71 Let A = {(a, 77, £) G A x X 2 : 77 < £} and assume 
A := V(A) G (0, 1). Since the coupling defined by L is monotone, we have 
l A (a t ,r] t ,£ t ) > l A (a ,?7 ,^o)- By Lemma E31 

z7 A := F(da, c?r/, d^|(a, 77, () 6 A) G X^ 

From 77 = Xu A + (1 — A)Fyic, we deduce zT^ £ X^. Since A is shift invariant 
in A x X 2 , V A and 174c lie in S. By extremality of V, we must have V A = V A c 
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which is impossible since these measures are supported on disjoint sets. □ 

Attractiveness assumption (A5) ensures that an initially ordered pair of 
coupled configurations remains ordered at later times. Assumptions (Al), 
(A4) induce a stronger property: pairs of opposite discrepancies between two 
coupled configurations eventually get killed, so that the two configurations 
become ordered. 

Proposition 3.2 Every V £ fl S is supported on O. 

Proof of proposition \3.2\ We follow the scheme used in [231 H El EE3 EES] 
for the non-disordered case, and only sketch the arguments needed for the 
disordered setting. 

Step 1. For x G Z, let f x {v,£) — (v( x ) ~~£,( x )) + - By translation invari- 
ance of u, the shift commutation property (fT9j) and (l2"Tj) . (|28|) . 



Using Assumptions (A4)-(A5), (JI}) and translation invariance of u, we obtain 



for x 7^ y with p(y — x) + p(x — y) > 0. Whenever one of the events in (1331) 
holds, we say there is a pair of opposite discrepancies at (x, y). 

Step 2. One proves by induction that, for all n € N, (1331) holds if x 7^ y 
with p* n (y — x) + p* n (x — y) > 0. The induction step is based on the fol- 
lowing idea. Assume (77, £) has a pair of opposite discrepancies at (x,y). 








On the other hand (see [HI [15] ) 




+ F((a,r/,0 : r)(x) < £(x), r](y) > £(y)) = 



(33) 
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Then one can find a finite path of coupled transitions (with rates uniformly 
bounded below thanks to (A4)-(A5) and (pQ)), leading to a coupled state with 
a pair of opposite discrepancies, either at (x, z) for some z with 
x) + p*^- 1 ) (x - z) > 0, or at (z, y) with p< n ^ (y — z) + p*^ 1 ) (z — y) > 0. 
This part of the argument is insensitive to the presence of disorder so long 
as a(x) is uniformly bounded below. 

Conclusion. By irreducibility assumption (Al), ( 133]) holds for all (x,y) G Z 2 
with x 7^ y. This implies F((9) = 1. □ 

We are now in a position to prove Proposition 13.11 

Proof of proposition \3.1\ We define 

1Z Q := |y v(Q) u (da, drf) : f G (Zl fl <S) e , i> has a-marginal Qj 

Let z/ G (X L Pl<S) e with a-marginal Q and p* := J r/^u^da, drj) G 72.^ for 
i G {1,2}. Assume p 1 < p 2 . Using Lemma \3.1\ (iii). Lemmas 13.21 and 13.4] 
and Proposition 13.21 we obtain v 1 <C z/ 2 , that is 023]). Existence ( 12~2]) of 
an asymptotic particle density can be obtained by a proof analogous to [25] 
Lemma 14], where the space-time ergodic theorem is applied to the joint 
disorder-particle process. Then, closedness of 1Z Q is established as in jU 
Proposition 3.1]. We end up proving the weak continuity statement given 
the rest of the proposition. Let p, p' G HP with p < p' . By ( J23l) and Lemma 
I3.1j there exists a coupling J7Q >p>p (da, dr] 1 d^) of u^' p (da, drj) and u Q ' p '(da : d£) 
supported on + . Thus, for x G Z 




(34) 



from which weak continuity follows by a coupling argument. 



□ 



Remark 3.3 Since 




(35) 



a Lipschitz constant V for G Q follows from /[23\) , (34\): 




(36) 
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4 Proof of hydrodynamics 



In this section, we prove the hydro dynamic limit following the strategy in- 
troduced in [3J H] and significantly strengthened in [5]. That is, we reduce 
general Cauchy data to step initial conditions (the so-called Riemann prob- 
lem) and use a constructive approach (as in [2J). Some technical details 
similar to [5] will be omitted. We shall rather focus on how to deal with the 
disorder, which is the substantive part of this paper. The measure Q being 
fixed once and for all by Theorem 12.11 we simply write u p , TZ, G. 

4.1 Riemann problem 

Let X, p G [0, K\ with A < p (for A > p replace infimum with supremum 
below), and 

R\, p {x, 0) = Al {;r< o} + pl{ x > } (37) 

The entropy solution to the conservation law (I13p with initial condition (|37|) . 
denoted by R\^ p (x,t), is given ([U Proposition 4.1]) by a variational formula, 
and satisfies 

/w 
R x , p (x,t)dx = t[G v / t (\,p) - G w /t{\p)\, with 

g v (X,p) := inf {G(r) - vr : r G [A, p] n 71} (38) 

for all v, w G R. Microscopic states with profile ( 1371) will be constructed 
using the following lemma, established in Subsection 14.31 below. 

Lemma 4.1 There exist random variables a and (r] p : p G TV) on a proba- 
bility space (^a^a^IPa) such that 

(a,^)~^ p , a~Q (39) 
Pa — as., p r] p is nondecreasing (40) 

Let V x ' p denote the distribution of (a,r] x ,ri p ), and u x ' p the conditional distri- 
bution of (a, r/ x , r] p ) given a. For (x , t ) £ Zx R + , the space-time shift 9 XO}to 
is defined for any uo G Q, for any (t, x, z, u) G M + x Z x Z x [0, 1], by 

(t, x, z, u) G XOtto co if and only if (t + t, x + x, z,u) G u) 



16 



By its definition and property ([70, the mapping introduced in (jlj) satisfies, 
for all s, t > 0, x G Z and (a,t),u)6AxXx £7: 

r} s (a,r] t (a,r},u),6o,tu) = r] t+s (a,r],uj) 

T x r] t (a,r],uj) = 7] t (T x a,T x r],9 xfi uj) 

We now introduce an extended shift 6' on f2' = A x X 2 x £7. If uj' = (a, t], £, u) 
denotes a generic element of Q', we set 

6'x,t w ' = ( T x a > T xVt(oi, r], u), r x r] t (a, £, w), 6*^) (41) 

It is important to note that this shift incorporates disorder. Let T : X 2 — > X 
be given by 

T(rj, 0(a) = 7?(x)l {;c< o} + e(^)l{x>o} (42) 
The main result of this subsection is 

Proposition 4.1 Set, for t > 0, 

ff(u')(dx) : = n N (r] t (a, T(rj, g),u))(dx) (43) 
For all t > 0, sq > and Xq G R, we have that, for Q-a.e. a G A, 
\im^ t (e[ NxoLNs y)(dx) = R x>p {.,t)dx, V x a p ® P-a.s. 

Proposition 14.11 will follow from a law of large numbers for currents. Let 
x_ = (x t , t > 0) be a Z-valued cadlag random path, with \x t — x t -\ < 1, 
independent of the Poisson measure u. We define the particle current seen 
by an observer travelling along this path by 

<Pt-(a,ijo,u) = (ft-' + (a,r] ,u) - (ft-~(a,r] ,u) + (pf(a,r}o,u) (44) 

where (p*'' (a,r)o,u) count the number of rightward/leftward crossings of x, 
due to particle jumps, and <f^ (a,r]o, u>) is the current due to the self-motion 
of the observer. We shall write tp\ in the particular case x t = [vt\ . Set 
0JV) := ^(a, Tfa, that for {v,w) G R 2 , ^ t (tu')([v,w]) = 
t(A^^)- 1 (^ / ^ ^ (a; , ) - <$f(oS)). By §E\), Proposition O is reduced to 

Proposition 4.2 For allt>0, a G R + , b G R and !)6l, 

lim (iVt)-V5r«( Vj.aiX) = &(A, p) ® P - a.s. (45) 



17 



To prove Proposition 14. 2\ we introduce a probability space Q + , whose generic 
element is denoted by u + , on which is defined a Poisson process (N t (oj + )) t > 
with intensity \v\ (v £ K). Denote by P + the associated probability. Set 

x?(u+) := (sgn(v)) [N aN+s (u+) - N aN (oj+)] (46) 
rj"(a,r] ,u,u + ) := T x N {uj+) r] s (a } ri ,uj) (47) 
5f(a,w + ) := T x N {bJ+) a (48) 

Thus (a^f ,rj^) s >o is a Feller process with generator 

U = L + S v , S v f(a, C) = \v\ [f(r sga[v) a, T sga(v) () - /(or, ()) 

for / local and a £ A, £ £ X. Since any translation invariant measure on 
A x X is stationary for the pure shift generator S v , we have II H5 = n«S. 
Define the time and space-time empirical measures (where e > 0) by 

m tN (uj',u + ) := (Nty 1 / <y(SAr(a^+),^(a,r(^o,w,w+))ds (49) 
m tN)£ (u', u + ) : = iZnt-e^eiVir 1 ^ r x m tN {u' (50) 

xeZ: |a;|<ejV 

Notice that there is a disorder component we cannot omit in the empirical 
measure, although ultimately we are only interested in the behavior of the 
77-component. Let -Ma,p denote the compact set of probability measures 
fi(da, drj) £ X L fl S such that /1 has a-marginal Q, and u x <C \i <C v p . By 
Proposition 13. 1[ 

M x, P = ^u(da, dr)) = J u r (da, dr))j(dr) : 7 £ V{[\, p) n K) J (51) 

The key ingredients for Proposition 14.21 are the following lemmas, proved in 
Subsection 14.31 below. 

Lemma 4.2 The function <f>t(a,r),£,u) is increasing in r\, decreasing in £. 

Lemma 4.3 With V x ' p ® P <S> P + -probability one, every subsequential limit 
as N -)> 00 o/m iA r j£ (6 l ' bAr j iaAr c<; / ,c<; + ) Zzes zn -Ma, p - 
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Proof of Proposition \4-2\ We will show that 

^M(Nt)- 1 <f> v tN od', bNi aN (u') > G v (X,p), ® P-a.s. (52) 

limsup (iV*)-Vrjv°0Uv;v(y) < &(A,p), P-a.s. (53) 

5"tep one: proo/ of ( tH| ). 

Setting ro aiV = cc7 aJ v(w') := T (r LbA rj?7 a7V (a, 77, w), r LbA rj?7 aA r(a, f, w)), we have 

(iVt) _ VtJV ° ^L&iVJ^JvC^O = (^) _ Vw(nWJ«' OlbN\,aNU) (54) 

Let, for every (a, (, w, u + ) GAxXxUxfl 1 and x N (u + ) given by 

V$(a,C, W ,a; + ) := |Z n [-eJV, ^] p 1 £ y^ +)+ >,C,a,) (55) 

y£Z: |y|<eV 

Note that lim A r_ 5 . 00 (Aft) _1 2^ r (c<; + ) = v, P + -a.s., and that for two paths y,,z 
(see dSD), 

IVwi^Vo,^) - <Pw(u,Vo,u)\ < K(\y tN - z tN \ + \y - z \) 

Hence the proof of (152]) reduces to that of the same inequality where we re- 
place (Nt)- 1 ^ o 9' [bN]AN (uj') by (M)- 1 ^^ 1 "^") ro «jv, ^wj.ajvw, w + ) and 
V x ' p <g> P by V x ' p <g) P <g> P + . By definitions dHJ), (jSJ) of flux and current, for 
any a G A, ( G X, 

tjv (a>C>w,w + ) := <^r (a,C,w)- 

/ ^ {j("f (a, (a, C, ^, - v (Vs( a i C, w, w + ))(l {tI>0} )} 

./o 

is a mean martingale under P eg) P + . Let 

B% := (ivtiznheiv^iv]!)- 1 £ M t ^ 
= (Art)- 1 ^(r L6JV jQ; ) K7 aiV) L6iV j iaJV a;, w + ) 

~y [i(a^) - v v( 1 {v>o})]m t N, e (0[bN\,aN U} ' ',u + )(da,dr]) (56) 
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where the last equality comes from (15 Op . (j55p . The exponential martingale 
associated with yields a Poissonian bound, uniform in (a,(), for the 

exponential moment of with respect to P <S> P + . Since zu a N is indepen- 
dent of (9[pN\,aN<^, w+ ) under I7 A ' P <g> P <g> P + , the bound is also valid under 
this measure, and Borel-Cantelli's lemma implies limjv-Kjo -R^v = 0. From 
( 151)]) . Lemma S3] and Corollary 13.14 (ii) imply f ]5"2"]) . as well as 

limsup(iVt)-Vijv ^j,aJv(w') < sup [G(r)-vr], V x > p ® P-a.s. (57) 

Step two: proof of [53) . Let r e [A, p] DTZ. We define z7 A ' r '^ as the distribution 
of (a,r] x ,r] r ,r]P). By and (J5ZD, 

lim (M)" ViW ° Vj,aiv(«> ^ ^ w) = G?(r) - w 
By Lemma [4.21 

$W ° e [bN\,aNW) < <PtN ° 0LWVJ,aAr(«> ^ ^ W ) 

The result follows by continuity of G and minimizing over r. □ 
4.2 Cauchy problem 

For two measures /i, z/ 6 A1 + (M) with compact support, we define 

A(/x, z/) := sup |z/((— oo, x]) — oo, x])| (58) 

which satisfies: //Pij For a sequence (/i„)„>o of measures with uniformly 
bounded support, //„—)■// vaguely is equivalent to lim n ^. 0O A(/i n , /i) = 0; 
(P2) the macroscopic stability property ([TJ [25] ) states that A is, with high 
probability, an "almost" nonincreasing function of two coupled particle sys- 
tems; (PS) correspondingly, there is A-stability for ( !T3|) . that is, A is a 
nonincreasing function along two entropy solutions ([5] Proposition 4.1, Hi), 

b)})- 

Proposition 4.3 Assume (t]q) is a sequence of configurations such that: 
(i) there exists C > such that for all N e N, rj^ is supported on Z R 
[-CN.CN]; (ii) 7T N (t]q) — > uo(.)dx as N — >• oo, where uo has compact 
support, is a.e. TZ-valued and has finite space variation. Let u(.,t) denote 
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the unique entropy solution to < f73j) with Cauchy datum Uq(.). Then, Q®P- 
a.s. as N — )■ oo, 

A"(f) := A(7r JV «(a,^, W )), M (.,t)rfx) 
converges uniformly to on [0,T] /or every T > 0. 

Theorem 12.11 follows for general initial data Uq by coupling and approxima- 
tion arguments (see [5J Section 4.2.2]). 

Proof of proposition \4 . 3\ By initial assumption f|T2|) . liniAr^oo A Ar (0) = 0. Let 
e > 0, and e 7 = e/(2V), for V given by §k>. Set * fc = he' for fc < « := [T/e'J , 
t K+ i = T. Since the number of steps is proportional to e, if we want to bound 
the total error, the main step is to prove 

limsup sup [A N (t k+1 )-A N (t k )]<38e, Q <g> P-a.s. (59) 

iV^oo fc=0,...,/C-l 

where 5 := 6(e) goes to as e goes to 0; the gaps between discrete times 
are filled by an estimate for the time modulus of continuity of A N (t) (see 
Lemma 4.5]). 

Proof of / f5Pj) . Since u(.,t k ) has locally finite variation, by [5j Lemma 4.2], 
for all e > we can find functions 

h 

Vk = J2 Tk > ll lWk,i+i) ( 60 ) 
1=0 

With -OO = £fc, < < • • • < < X k ,l h +1 = +°°> r M e ^> r fc,o = ?"fc,Z fe = 

0, such that x k ,i — x k ,i-i > e, and 

A(u(.,t k )dx,v k dx) < 8e (61) 

For t k <t < tfc+i, we denote by v k (., t) the entropy solution to f|T3|) at time t 
with Cauchy datum v k (.). The configuration £ N,k defined on (Qa. ® ^, <8> 
J 7 , P A <S> P) (see Lemma EQ by 

f • fe (u A) u)(x) := i]jvt t (a(uA),1 rM N,w)(i), if L^fc,iJ < & < L^fc,i+iJ 
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is a microscopic version of v k (-), since by Proposition 14. II with A = p = r k 

,N/<-N,k 



lim /(f lfc (w A ,w))(dx)=^(.)(lx, P A ® P-a.s. 

N— >oo 



(62) 



We denote by ^ t ' (wa,w) = 77 t (a(u;A), £^ ,fc (^A, w), 0o,m fc ^) evolution start- 
ing from £ N ' k . By triangle inequality, 

A N (t k+1 )-A N (t h ) < A[/(i +i ),/0]-A w (4) (63) 



A 



+ A(v h (.,s')dx,u(.,t k+1 )dx) 



(64) 
(65) 



To conclude, we rely on Properties (Pl)-(PS) of A: Since e' = e/(2V), finite 
propagation property for (fT3"]l and for the particle system (see [H Proposition 
4.1, Hi), a) and Lemma 4.3]) and Proposition 14.11 imply 

lim 7r N (£ N ',(u) A , u)) — v k (.,£') dx, ® P-a.s. 

N— >oo 

Hence, the term (IMj) converges a.s. to as N — > oo. By A-stability for (fl"3"j) . 
the term (165]) is bounded by A(v k (-)dx, «(., t k )dx) < 5e. We now consider 
the term f )63|) . By macroscopic stability (j25l Theorem 2, Equation (4) and 
Remark 1]), outside probability e~ CNSe , 



A 



* N (v% tk J, * N (& < A 7r"(£**)] + fe (66) 



Thus the event ( 166]) holds a.s. for N large enough. By triangle inequality, 

aK(%\),^(H-^W 

< A («(., * fc )dz, u fc (.)da;) + A [v fc (.)dx, 7r"(f 

for which flBTj) . ( 162|) yield as iV -> oo an upper bound 25e, hence 3fe for the 
term (J63J. □ 

4.3 Proofs of lemmas 



Proof of lemma J^.l. Let 1Z<i be a countable dense subset of 1Z that con- 



tains all the isolated points of 71. We denote by TZt, resp. TZ d , the set 
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of p G [0, K] that lie in the closure of [0,p) fl IZj, resp. (p,K] R 7^. Be- 
cause 7Z is closed, we have 7Z = TZd U 7£j U 7Z d . By (I2H|1 there exists a 
subset A' of A with Q-probability 1, such that u p < u p for all a G A' and 
p, p' G 7^- By [19, Theorem 6], for every a G A', there exists a family of 
random variables [r\ p a : p G on a probability space (O a , J a ,P a ), such 
that dHH]) (HOD hold for p G 7t d . Let f2 A = {(a,u a ) : a G A', w Q G 
J-"a be the a-field generated by mappings (a,co a ) h-> p p (a:,c<; Q ) := ^(a; a ) for 
p G 7^d, and PA(da,dco a ) = Q(da) g) P a (cL; Q ,). Now consider p E 7Z\ TZd- 
Since p^ is a nondecreasing function of r, for every a G A' and u a G fi a , 
p p+ (a,w Q ) := lim r ^ Pjr<Pire ^ d r/^(w a ) exists if p G 72^, and p p ~(a,u; Q ) := 
lim r ^ Pir>p>re ^ d ^(w a ) exists if p G 7Z d . We set r) p (a,u a ) = r] p+ (a,u a ) if 
p G T^J, p p (a,c<; a ) = ^"(a,^) otherwise. Suppose for instance p G 7£j. 
Since p p+ is a PA-a.s. limit of p r as r — )■ p, r < p, r G 7^, it is a limit in 
distribution. Weak continuity of v p then implies ( 139]) . Property fj40|) on TZ 
follows from the property on TZd and definitions of r] p± . □ 

To prove Lemma 14.3} we need the following uniform upper bound (proved in 
[3 Lemma 3.4]). 

Lemma 4.4 Let ~P V V denote the law of a Markov process (5.,£.) with gener- 
ator L v and initial distribution v. For e > 0, let 

7r ti£ := |Zn \-etM\~ 1 E W ^.S-.r-S)* 8 ( 67 ) 

xeZn[-£t,et] 

Then, there exists a functional 7) v which is nonnegative, l.s.c, and satisfies 
D^ 1 (0) = Xu», such that, for every closed subset F ofV(A x X), 

limsup/TMog sup P v u U t ,e(0 e F) < - inf V v (p) (68) 

Proof of Lemma \4-3\ We give a brief sketch of the arguments (details are 
similar to [51 Lemma 3.3]). Spatial averaging in (157)) implies that any subse- 
quential limits p lies in 5. Lemma I4~4l and Borel-Cantelli's Lemma imply that 
p lies in X^v (uniformity in ( 168)) is important because ^-shifts make the initial 
distribution of the process unknown). Finally, the inequality v x <C p <C v p 
is obtained by coupling the initial distribution with r/ x and rj p , using 
tiveness and space-time ergodicity for the equilibrium processes. □ 
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Proof of lemma \4~J^ Assume for instance rj < rf . Let 7 := T(j},£) and 
i ■= T(if,ti), It = r] t (a,-f,uj) and j' t = r] t (a,i ',u). By (JSJ), j t < j' t for 
all t > 0. By definition of the current, (p v t (a, 7/, £, lo) — (p^(a,r],^,uj) = 

j: x> m(x) - it(x)} > 0. ' □ 

5 Other models 

For the proof of Theorem 12.11 we have not used the particular form of L a in 
(T5]), but the following properties. 

1) The set of environments is a probability space (A,J-a ; Q) ; where A is 
a compact metric space and JFa its Borel cr-field. On A we have a group 
of space shifts (t x : x G Z), with respect to which Q is ergodic. For each 
a G A, L a is the generator of a Feller process on X that satisfies (fT9|) . The 
latter should be viewed as the assumption on "how the disorder enters the 
dynamics". It is equivalent to L satisfying (1T51) . that is being a translation- 
invariant generator on A x X. 

2) For L a we can define a graphical construction |5| ona space-time Poisson 
space P) such that L a coincides with fTiTJj) . for some mapping T a,z,v 
satisfying the shift commutation and strong attractiveness properties ([?]) and 
([8]). The existence of this graphical construction for the infinite- volume sys- 
tem follows from assumption (A2), which controls the rate of faraway jumps. 
This assumption is also responsible for the finite propagation property of 
discrepancies in the particle system, and its macroscopic counterpart, the 
Lipschitz continuity of the flux function (see ( TTBTl . Remarks 13.21 and 13. 3p . 

3) Irreducibility and non-degeneracy assumptions (Al), (A4) (combined with 
attractiveness assumption (A5)) imply Proposition 13.21 

In the sequel we consider other models satisfying 1) and 2), for which ap- 
propriate assumptions replacing (A1)-(A5) imply existence of a graphical 
construction, and Proposition 13.21 as in 3). In these examples, the transition 
defined by 7 -a ' 2,t ' in (JSj) is a particle jump, that is of the form T a < z < v ri = 



24 



^(o.*.*),v(o.*.«). It follows that (HDD yields (in replacement of P» 
where 

c Q (x, y, rj) = m (i v e V : T a ' z <^ = 77**}) (70) 

and the shift-commutation property (j7|) implies 

c a (x,y,T]) = c Txa (0,y- x,r x rj) (71) 

which, for ( )69|) . is equivalent to ( TT9|) . Microscopic fluxes ( JT4j) and (125"]) more 
generally write 

j + (a,r)) = c a(v(y),v(y + z)) 

y,z£Z: y<0<y+z 

j~(a,rj) = c a{v(y),v(y + z)) 

y.z(r.%: y+z<0<y 

J(a,v) = ^2zc a (0,z,r]) (72) 

zez 

5.1 Generalized misanthropes' process 

Let c G (0,1), and p(.) (resp. -P(-)), be a probability distribution on Z 
satisfying assumption (^4ij (resp. (A2j). Define A to be the set of functions 
B : 1? x {0, . . . , K} 2 M+ such that for all (x, z) e Z 2 , B(ac, 2, ., .) satisfies 
assumptions (A3)-(A5) and 

50,2,1,^-1) > cp(» (73) 
z,#,0) < c _1 P(z) (74) 

The shift operator r y on A is defined by (r y B)(x, z, n, m) = B(x + y, z, n, m). 
We generalize fl2]) by setting 

L a f(v) = B & y - x > ^1 V(y)) [f (V x ' y ) - f{v)} (75) 

where we assume that the distribution Q of B(., .,.,.) is ergodic with respect 
to the above spatial shift (we kept the notation L a to be consistent with 
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the rest of the paper, but we should have written L B ). Assumption ( !73"j) 
replaces (Al) and implies Proposition 13.21 Assumption (J7I|) replaces (A2) 
and implies existence of the infinite volume dynamics given by the following 
graphical construction. For v = (z,u), set m(dv) = c~ 1 P(dz)X[ ^](du) in fl3]), 
and replace with 

+z B{x 1 z 1 r ] {x) 1 r 1 {x + z)) 

T a ^ v r 1 ={ 71 U< c - l P{z) (76) 

7] otherwise 



Here the microscopic flux (1721) writes 

and the Lipschitz constant V = 2c _1 ^2 z€li \z\P(z) for follows as in 
fl36|) from fl34|) - fl35l) . The basic model (T5]) is recovered with B(x, z,n,m) = 
a(x)p(z)b(n,m), for p(.) a probability distribution on Z satisfying (A1)-(A2), 
ct(.) an ergodic (c, l/c)-valued random field, and .) a function satisfying 
(A3)-(A5). In this case (f73]) - (f71]) hold with P(.) = Here are two other 
examples. 

Example 1.1. This is the bond-disorder version of ([2]): we have B(x, z, n, m) = 
a(x,x + z)b(n,m), where a = (a(x,y) : x, y G Z) is a positive random 
field on Z 2 , bounded away from 0, ergodic with respect to the space shift 
r z a — a(. + z, . + z). Sufficient assumptions replacing (Al) and (A2) are 

cp(y — x) < a(x, y) < c~ 1 P(y — x) (77) 

for some constant c > 0, and probability distributions p{.) and P(.) on Z, 
respectively satisfying (Al) and (A2). 

Example 1.2. This is a model that switches between two rate functions ac- 
cording to the environment: we have B(x, z, n, m) = p(z)[{l— a(x))b (n, m) + 
a(x)bi(n,m)}, where (a(x), x G Z) is an ergodic {0, l}-valued field, p(.) sat- 
isfies assumption (Al), and b , b\ assumptions (A3)-(A5). 



5.2 Generalized /c-step A'-exclusion process 

We first recall the definition of the fc-step exclusion process, introduced in 
[TT] . Let K — 1, k G N, and p(.) be a jump kernel on Z satisfying assump- 
tions (A1)-(A2). A particle at x performs a random walk with kernel p(.) 
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and jumps to the first vacant site it finds along this walk, unless it returns 
to x or does not find an empty site within k steps, in which case it stays at 
x. 



To generalize this, let K > 1, k > 1, c G (0,1), and T> denote the set of 
functions (3 = . . . , (3 k ) from Z k to (0, l] fc such that 

PK) e [c,l] (78) 
P%) > /3 l+1 (.),V*G{l,...,fc-l} (79) 

In the sequel, an element of Z fc is denoted by z = (z±, . . . ,Zk). Let q be a prob- 
ability distribution on Z fc , and (3 G V. We define the (q, (3)-k step .fT-exclusion 
process as follows. A particle at x (if some) picks a g-distributed random vec- 
tor Z_ = (Zi, . . . , Zk), and jumps to the first site x + Z{ (i G {1, . . . , k}) with 
strictly less than K particles along the path (x + Z\, . . . , x + Zk), if such a 
site exists, with rate (3 l (Z_). Otherwise, it stays at x. The /c-step exclusion 
process corresponds to the particular case where K — 1, q is the distribution 
(hereafter denoted by q^wip)) of the first k steps of a random walk with 
kernel p(.) absorbed at 0, and (3 l (z) = 1. Outside the fact that K can take 
values > 1, our model extends fc-step exclusion in different directions: 

(1) The random path followed by the particle need not be a Markov process. 

(2) The distribution q is not necessarily supported on paths absorbed at 0. 

(3) Different rates can be assigned to jumps according to the number of steps, 
and the collection of these rates may depend on the path realization. 



Next, disorder is introduced: the environment is a field a = ((q x ,(3 x ) '■ x G 
Z) G A := (V(Z k ) x T>)' L . For a given realization of the environment, the 
distribution of the path Z_ picked by a particle at x is q x , and the rate at 
which it jumps to x + Zi is (3 X {Z_). The corresponding generator is given by 
( |69l) with c a = J2i=i c «' where (with the convention that an empty product 
is equal to 1) 



c l a (x,y,ri) = 1 



{ V { X )>0}i-{r,{y)<K} 



i-1 



P x (z)l{x+z z =y} Y\_ 1 {V^ 



+ Zj )=K} 



dq x (z) 



The distribution Q of the environment on A is assumed ergodic with respect 
to the space shift r y , where T y a = {{q x+y , (3 x+y ) : x G Z). 
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For the existence of the process and graphical construction below, and for 
Proposition I3.2[ sufficient assumptions to replace (Al)-(A2) are 



x&L 



> cp(.) 



sup supg*,(.) < c P(.) 

i=l,...,k xdTL 



(80) 
(81) 



for some constant c > 0, where q x denotes the i-th marginal of q x , and p(-), 
resp. -P(.)> are probability distributions satisfying (Al), resp. (A2). To write 
the microscopic flux and define a graphical construction, we introduce the 
following notation: for (x,z,T]) G Z x Z fc x X, G V and u G [0,1], 



N(x, z, rf) 
Y(x,z,rj) 



inf {i G {1, . . . , k} : r/ (x + z«) < i^} withinf I 

X + Z N ( Xj zm) if N{x,z_,7]) < +oo 
if N(x,z, r]) = +oo 



+oo 



v x,Y(x,z,n) if ^ > and « < /^(^ (z) 
77 otherwise 



(where the definition of (3 +OQ (z) has no importance). With these notations, 
we have 



l{ ?? (x)>0}E (?0 



N(x,Z,7)) . 



-{r,(0)>0}E go 







-{Y(x,Z, V )=y} 



c a (x,y,v) 

where expectation is with respect to Z. Since 

1=1 i=l 

§M-§B y iel d f o r G Q the Lipschitz constant V = 2k 2 C - 1 Y. ze z 



(82) 
(83) 



)-i{Zi] 

\z\P(z). 



Let V = [0,1] x [0,1], m = A[o,i] <S> A[o,i]. For each probability distribu- 
tion q on Z fc , there exists a mapping F g : [0, 1] — > Z fc such that F q (Vi) has 
distribution q if Vi is uniformly distributed on [0, 1]. Then the transformation 
T in ()5]) is defined by (with v = (vi, v 2 ) and a = ((q x , (3 X ) : x G Z)) 



<-ra.,x,v _ q-^,F qx (vi),l3 x (F qx (v 1 )),V2 
I '/ — '0 '/ 

Strong attractiveness of our process will follow from 



(84) 
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Lemma 5.1 For every (x,z,u) G Z x Z fc x [0,1], j^'-'P' u i s an increasing 
mapping from X to X. 

Proof of lemma [3771 Let (77, £) G X 2 with 77 < £. To prove that 7^*% < 
7^'-^'"C) since 77 and £ can only possibly change at sites x, y := 77) 
and t/' := y(:r, 2, £), it is sufficient to verify the inequality at these sites. 

If £(x) = 0, then by ( |84l) . 7/ and £ are both unchanged by j^ &l3 ' u _ if 
r,(x) = < £(*), then V'" u XUf) > £(t/) > V (y') = X^iv')- 

Now assume r](x) > 0. Then 77 < £ implies N(x,z,rj) < N(x,z,^). If 
N(x,z,r)) = +00, 7/ and £ are unchanged. If N(x,z,r)) < N(x,z,£) = +00, 
then J^'-' 13 ^^ = rf^y an d £(7/) = i^. Thus, Tq u r](x) = r)(x) — 1 < £(2;) = 
7^^'"£(x) and V^Cl/) = Civ) = K > '77i' : '""'/(//!• If iV(x,z,r/) = 
< th en /^'^ = /^(^ =: /3. If u > p both 7/ and 

£ are unchanged. Otherwise Tq u 7] = ^ x >v an d j'^ , - ,l3 ' u ^ = £ a; >2/ ) whence 
the conclusion. Finally, assume N(x,z,r)) < N(x,z,£) < +00, hence (3 : = 
^au) > 0N(x,z,O = . p by ([Tgj) and r)( y ) <£( y ) = K. If 7j > /3 , 77 and £ are 

unchanged. If 77 < /?', then Tq u t](y) = r)(y) + 1 < £(7/) = Tq u ^(y) = K 
and ^'^'"£(7/') = £(t/') + 1 > K < « < A then V'^^ = 

77(x)-l < 7^' u £(*) and T x ^ u V (y) = V (y)+1 < T?'^(y) = £(y) = 
□ 

We now describe a few examples. 

Example 2.1. Let K = 1, (a x : £ G Z) be an ergodic [c, l/c]-valued ran- 
dom field, and r(.) be a probability measure on Z satisfying (A1)-(A2). A 
disordered version of the /c-step exclusion process with jump kernel r is ob- 
tained by multiplying the rate of any jump starting from x by a x . This 
means that the random field (q x ,(3 x ) xe z is defined by q x = qR W (r), and 
AcU) = (a x , ...,a x ) for every z G Z fc . 

Example 2.2. Let (7^, i^ez be an ergodic [c, l] 2fc -valued random field, where 
lx = (7™, 1 < 77 < /c) and l x = (7.™, 1 < 77 < k). The random field (q x , f3 x ) xe z 
is defined by 

= -0(1,2,...,*) + -2,... ,-fc) 
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#(1, 2, . . . , k) = 27;, #(-1, -2, • • • , -*) = 24 

Hence the rates are disordered but not the distribution of the random path 
followed by particles: the stationary random field (q x ) x &z is deterministic and 
uniform. Here, the jump rate and microscopic flux fl82l) - fl83l have a fairly 
explicit form: 

c a (x,y,r]) = '~?y- x l Mx)>0} l {ri{y)<K} Yl l{r,(z)=K} if V > x (85) 

z=x+l 
x-1 

c a (x,y,T]) = L x ~ v l {v{x)>0] l {n{y)<K] Yl 1 Mz)=K} if y < x (86) 

z=y+l 
k n—1 

J{ot,rj) = 77(0)^/17^(1 -rj{n)) J[vU) 

n=l j=l 
k n—1 

- v(0)J2<(l-v(-n))HvH) (87) 

71=1 j=l 

Example 2.3. Set q x = qKw(r x ), for (r x ) x£ z an ergodic random field with 
values in the probability measures on Z satisfying (Al)-(A2). The simplest 
case is nearest-neighbor jumps, that is, r x = p x 5i + (1 — p x )5-i, where, 
for some c G (0,1), (p x ) xt =z is an ergodic [c, l/c]-random field. Due to the 
nearest-neighbor assumption, a particle starting from x can only jump to 
y > x (resp. y < x) if y is not full and all sites between x and y (resp. y and 
sc) are full. Hence, the jump rate (J8~2"j) is identical (see example below) to the 
one obtained by taking in 



L(fc-n)/2j L(*-")/2j 

l U x= E P" + '(l-^n(n + /,/), 4= E (l-^) n+ ^n(n + M) 
/=0 i=0 

for u G {1 fc}, where C n (i,j), for i, j G Z + and i + j > 0, is the number 

of paths (zo = 0, . . . , Zi+j) such that < z m < n for m = 1, . . . , i + j — 1, 
l^m+i — z m \ = 1 for m = 1, . . . , i + j, and Card{m G {1, . . . , i + j} : z m — 
z m _i = 1} = i. With this choice of 7" and the microscopic flux is given 
by (1H7|) . For instance if k = 5, we obtain, for n G {1, . . . , k}\ 

n-1 

c a (x, x + n,r}) = p™l {v{x)>{i} l {ri{x+n)<K} Yl l{ v ( x +j)=K} if n ^ 3 (88) 

j'=i 
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c a (x,x + 3,?7) 



pl[l + p x (l - p x )} x 

T)(x)r)(x + l)r)(x + 2)(1 — i](x + 3)) 



(89) 



n-l 



c a (x, x — n,r)) 



(1 — Px) n ^-{ n {x)>0}^-{r){x-n)<K} II l{r ? (x-j)=_ff } 



if n ^ 3 



c Q (z, x - 3,?7) 



(1 -p x ) 3 [l +p x (l -p x )] x 



rj(x)rj(x — l)r)(x — 2)(1 — 77(21 — 3)) 



Indeed, for n > and n^3, the only path from x to x + n that reaches x + n 
in at most k steps before returning to is x, x + 1, . . . , x + n. For n = 3, the 
additional path x, x + 1, x + 2, x + 1, x + 2, x + 3 yields the factor ^(1 — p x ) 
in (I8"9"j) . For n < 0, we change p x to 1 — p x . 

Note that in this process a given particle does not follow a random walk 
in random environment (RWRE) before it finds a non full site, but a ho- 
mogeneous random walk depending (randomly) on its initial location. For 
instance, in a 3-step process, a particle initially at x G Z will follow the path 
x, x + l,x+ 2,x+ 1 with probability — 

Example 2.4- The same random field (p x ) x( zz gives a different model if, at 
each transition, the selected particle follows a RWRE (X n ) n > with transition 
probabilities 

P(X n+1 = x + l\X n = x)= Px , P(X n+1 = x -l\X n = x) = l-p x (90) 

That is, we let q x be the distribution of (Xf — x, . . . , X£ — x), for (X%, 1 < 
n < k) a length Markov chain starting at x with transition probabilities 
( 190]) . There, unlike in Example 2.3 above, a particle initially atieZ follows 
the path x, x + 1, x + 2, x + 1 with probability p x p x+ \(l—p x+ 2)- The generator 
of this process is also identical to that of example 2.2, with 7" and of the 



7x — 7 n G°y '■ x <y < x + n), i n x = i n {p y : x - n < y < x) 
for some polynomial functions 7", t n : [0, 1]™ — > [0, +00), where n G {1, . . . , k}. 



form 
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5.3 i^-exclusion process with speed change and traffic 
flow model 



Let /C := {— k, ...,k} \ {0}, and a = ((v(x),j3*) : x G Z) be an ergodic 
[0, +oo) 2fc x (0, +oo)-valued field, where v(x) = (v z (x) : z G JC). We define 
the following dynamics. Set 



In configuration rj, if Z(a,x,rj) > 0, a particle at x picks a site y at random 
in Q(x,rj) with probability Z(a, x, rj)~ x v y ^ x {x), and jumps to this site at 
rate [3 X . If Z(a,x,rj) = 0, nothing happens. For instance, if v z (x) = 1, 
the particle chooses a vacant site uniformly. The corresponding generator is 
given by f )69|) . with 

C a (x,y,T)) = l{r 1 (x)>0}i{Z(a,x, V )>0}'ie(x,r ] )(y)Z(a,X,T]y 1 V y ^ x (x) 

Hence, the microscopic flux ( |25l) writes 



This process can be compared with a bond-disordered K- exclusion process in 
which a particle at x jumps to y with rate a(x, y) = v y - x (x). The difference is 
that in the latter, the particle could pick an occupied location, in which case 
the jump is suppressed. In the former, the particle first eliminates occupied 
sites and picks an empty site whenever there is at least one. This results in a 
speed change ^-exclusion process, that is the jump rate from x to y has the 
form c x ^ y (ri)l^ x ) >0 jl^y) <K y . To illustrate this, consider a nearest-neighbor 
example: we take K — 1, k — 1, Vi(x) = p(x) G [0, 1], = 1 — p(x). If 

sites x — 1 and x + 1 are free, in both processes the particle at x moves with 
rate (3 X to a site picked in {x — 1, x + 1} with probabilities p(x) and 1 — p(x). 
Now assume x + 1 is free and x — 1 occupied. If p(x) = 0, nothing happens 
in either process. If p(x) > 0, at rate the particle at x moves to x + 1 
in the speed change process, while in the bond-disordered process it moves 
to x + 1 with probability p(x) and attempts in vain to jump to x — 1 with 
probability 1 — p(x). 



e(x,r)) 
Z(a, x, rj) 




{y G Z : y-x G /C, 7?(y) < K} 



2£6(l,)j) 



77) = /3 1 l {j)(0) >o}^(a, 0, r]) 1 ^ 2 (0)l {r)(z)< ^ } 



32 



Assume K — 1, and consider the totally asymmetric case, where v z (x) = 
for z < 0. Recalling that the totally asymmetric exclusion process is a clas- 
sical simplified model of single-lane traffic-flow (without overtaking) where 
particles represent cars, the above model can be viewed as a traffic-flow model 
with maximum overtaking distance k. This is true also for Example 2.2 in 
Subsection 15.21 in the totally asymmetric setting i l x — 0, 1 < i < k. However 
in the latter model, an overtaking car has only one choice for its new position. 

Though it is not clear from this formulation, we can rephrase this dynamics 
as a 2/c-step model, which is thus strongly attractive by Lemma 15.11 To 
this end we take a random field of the form (3 X = . . . ,(3l), and define 
q x := q(v(x)), where q(v z : z G /C) is the distribution of a random self- 
avoiding path (Z\, . . . , Z 2 k) in K- such that 

nzx = y) = (9i) 

P(Z i = y\Z 1 ,...,Z i _. 1 ) = — ^ for2<«<2A; (92) 

For this model, assumption (1811) is always satisfied, while (18 Op reduces to 
the existence of a constant c > and a probability distribution p(.) on Z 
satisfying assumption (Al), such that 

inf v.(x) > cp(.) 
The link between the two models comes from 

Lemma 5.2 Assume (Z\, . . . , Z 2 k) ~ q(v z '■ z G fC). Let 6 be a nonempty 
subset of {z G K : v z ^ 0}, r := inf{z G {1, . . . , 2k} : Z { G 6} ; and Y = Z T . 
Then 

P(K = y) = l e (y); ''" 



Proof of lemma [Jjgj For all t > 2, let Ot_i be the set of self-avoiding paths 
(zi, . . . , of size t - 1 on K \ 9. For y G 6, by jSTJ-®, 



2 A- 



P(Y = y) = ^P(Z t = y,r = t) 



t=i 
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y) + 



2k 



E E 



zi, ■ ■ ■ , Z t -\ — z t -i) 



V, 



V 



t=2 ( Z1 ,...,z t -i)ee t -i 




Zl,...,Zt-l} z 



Cv. 



y 



where C is independent of y G 0, whence the result. 



□ 



Acknowledgments: K.R. was supported by NSF grant DMS 0104278. We 
thank BCM at TIMC - IMAG, Universites de Rouen, Paris Descartes and 
Clermont 2, and SUNY College at New Paltz, for hospitality. 

References 

[1] Andjel, E.D. Invariant measures for the zero range process. Ann. Probab. 
10 (1982), no. 3, 525-547. 

[2] Andjel, E.D., Vares, M.E. Hydrodynamic equations for attractive parti- 
cle systems on Z. J. Stat. Phys. 47 (1987), no. 1/2, 265-588. Correction 
to : "Hydrodynamic equations for attractive particle systems on Z" . J. 
Stat. Phys. 113 (2003), no. 1-2, 379-380. 

[3] Bahadoran, C, Guiol, H., Ravishankar, K., Saada, E. A constructive 
approach to Euler hydrodynamics for attractive particle systems. Ap- 
plication to /c-step exclusion. Stoch. Process. Appl. 99 (2002), no. 1, 



[4] Bahadoran, C, Guiol, H., Ravishankar, K., Saada, E. Euler hydrody- 
namics of one-dimensional attractive particle systems. Ann. Probab. 34 
(2006), no. 4, 1339-1369. 

[5] Bahadoran, C, Guiol, H., Ravishankar, K., Saada, E. Strong hydrody- 
namic limit for attractive particle systems on Z. Elect. J. Probab. 15 
(2010), no. 1, 1-43. 

[6] Benjamini, I., Ferrari, P.A., Landim, C. Asymmetric processes with ran- 
dom rates. Stoch. Process. Appl. 61 (1996), no. 2, 181-204. 



1-30. 



34 



[7] Bramson, M., Mountford, T. Stationary blocking measures for one- 
dimensional nonzero mean exclusion processes. Ann. Probab. 30 (2002), 
no. 3, 1082-1130. 

[8] Cocozza-Thivent, C. Processus des misanthropes. Z. Wahrsch. Verw. 
Gebiete 70 (1985), no. 4, 509-523. 

[9] Dai Pra, P., Louis, P.Y., Minelli, I. Realizable monotonicity for 
continuous-time Markov processes. Stock. Process. Appl. 120 (2010), 
no. 6, 959-998. 

[10] Evans, M.R. Bose-Einstein condensation in disordered exclusion models 
and relation to traffic flow. Europhys. Lett. 36 (1996), no. 1, 13-18. 

[11] Faggionato, A. Bulk diffusion of ID exclusion process with bond disor- 
der. Markov Process. Related Fields, 13 (2007), no. 3, 519-542. 

[12] Faggionato, A., Martinelli, F. Hydrodynamic limit of a disordered lattice 
gas. Probab. Theory Related Fields, 127 (2003), no. 4, 535-608. 

[13] Fill, J. A., Machida, M. Stochastic monotonicity and realizable mono- 
tonicity. Ann. Probab. 29 (2001), no. 2, 938-978. 

[14] Fritz, J. Hydrodynamics in a symmetric random medium. Commun. 
Math. Phys. 125 (1989), no. 6, 13-25. 

[15] Gobron, T., Saada, E. Couplings, attractiveness and hydrodynamics for 
conservative particle systems. Ann. Inst. H. Poincare Probab. Statist. 
46 (2010), no. 4, 1132-1177. 

[16] Gongalves, P., Jara, M. Scaling limits for gradient systems in random 
environment. J. Stat. Phys. 131 (2008), no. 4, 691-716. 

[17] Guiol, H. Some properties of /c-step exclusion processes. J. Stat. Phys. 
94 (1999), no. 3-4, 495-511. 

[18] Jara, M. Hydrodynamic Limit of the Exclusion Process in Inhomoge- 
neous Media. In: M. M. Peixoto; A. A. Pinto; D. A. Rand. (Eds.). 
Dynamics, Games and Science II. DYNA 2008, in Honor of Mauricio 
Peixoto and David Rand, University of Minho, Braga, Portugal, Septem- 
ber 8-12, 2008, Springer Proceedings in Mathematics, 2011, 2, 449-465. 



35 



[19] Kamae, T., Krengel, U. Stochastic partial ordering. Ann. Probab. 6 
(1978), no. 6, 1044-1049. 

[20] Kipnis, C, Landim, C. Scaling limits of interacting particle systems. 
Grundlehren der Mathematischen Wissenschaften [Fundamental Princi- 
ples of Mathematical Sciences], 320. Springer- Verlag, Berlin, 1999. 

[21] Koukkous, A. Hydrodynamic behavior of symmetric zero-range pro- 
cesses with random rates. Stochastic Process. Appl. 84 (1999), no. 2, 
297-312. 

[22] Krug, J., Seppalainen, T. Hydrodynamics and Platoon formation for a 
totally asymmetric exclusion model with particlewise disorder. J. Stat. 
Phys. 95 (1999), no. 3-4, 525-567. 

[23] Liggett, T.M. Coupling the simple exclusion process. Ann. Probab. 4 
(1976), no. 3, 339-356. 

[24] Liggett, T.M. Interacting particle systems. Classics in Mathematics 
(Reprint of first edition), Springer- Verlag, New York, 2005. 

[25] Mountford, T.S., Ravishankar, K., Saada, E. Macroscopic stability for 
nonfinite range kernels. Braz. J. Probab. Stat. 24 (2010), no. 2, 337-360. 

[26] Nagy, K. Symmetric random walk in random environment in one dimen- 
sion. Period. Math. Hungar. 45 (2002), no. 1-2, 101-120. 

[27] Quastel, J. Bulk diffusion in a system with site disorder. Ann. Probab. 
34 (2006), no. 5, 1990-2036. 

[28] Rezakhanlou, F. Hydrodynamic limit for attractive particle systems on 
Z d . Comm. Math. Phys. 140 (1991), no. 3, 417-448. 

[29] Seppalainen, T. Existence of hydrodynamics for the totally asymmetric 
simple i^-exclusion process. Ann. Probab. 27 (1999), no. 1, 361-415. 

[30] Serre, D. Systemes de lois de conservation. Diderot Ed., 1996. 

[31] Strassen, V. The existence of probability measures with given marginals. 
Ann. Math. Statist. 36 (1965), no. 2, 423-439. 



36 



